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Abstract 

We present X-ray observations of the northern outskirts of the relaxed galaxy 
cluster A1413 with Suzaku, whose XIS instrument has the low intrinsic background 
needed to make measurements of these low surface brightness regions. We excise 15 
point sources superimposed on the image above a flux of 1 x 10~^^ erg cm~^ s~^ (2- 
10 keV) using XMM- Newton and Suzaku images of the cluster. We quantify all known 
systematic errors as part of our analysis, and show our statistical errors encompasses 
them for the most part. Our results extend previous measurements with Chandra 
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and XMM-Newton, and show a significant temperature drop to about 3 keV at the 
virial radius, r2oo- Our entropy profile in the outer region (> 0.5 r2oo) joins smoothly 
onto that of XMM-Newton, and shows a flatter slope compared with simple models, 
similar to a few other clusters observed at the virial radius. The integrated mass of 
the cluster at the virial radius is approximately 7.5 x and varies by about 

30% depending on the particular method used to measure it. 

Key words: galaxies: clusters: individual (Abell 1413) — X-rays: galaxies: 
clusters — X-rays: ICM 



1. Introduction 

X-ray observations of intracluster medium (ICM) primarily give density and temper- 
ature. Density may be deduced relatively straightforwardly from cluster surface brightness 
because the ICM is optically thin and the emission coefficient over most observed bandpasses 
is nearly independent of temperature. There is good agreement on the ICM density over the 
observed radial range among different observers. In contrast, cluster temperatures have not 
been measured much beyond about half of the virial radius and, until recently, the shape of the 
temperature radial profile was a matter of heated debate even to that radius. Now independent 
measurements using four different observatories are consistent with a factor of ~ 2 decline of 
the projected temperature from the center to half the virial radius, at least for relaxed clusters 
(Markevitch et al. 1998; De Grandi & Molendi 2002; Vikhlinin et al. 2005; Piffaretti et al. 
2005; Pratt et al. 2007). 

In addition to the fundamental density and temperature variables, it is possible to derive 
additional thermodynamic variables from them, such as pressure and entropy. These derived 
variables are very useful when trying to understand the gravitational and non gravitational 
processes that were operative during the formation and evolution of the cluster. With the 
assumption of hydrostatic equilibrium, the cluster's total mass can also be derived from the 
ICM temperature and the radial derivatives of temperature and density. While this assumption 
is not valid for many clusters. X-ray observations are one of the few that can be used to measure 
such masses. 

The cold dark matter (CDM) paradigm combined with numerical simulations predict 
that the structure of clusters should exhibit self-similar scaling. That is their properties should 
be the same when scaled appropriately by redshift and the virial radius. This expected behavior 
occurs because clusters form from scale-free density perturbations and their evolution is mainly 
set by scale-free gravity, and both of these result because cluster masses are mainly CDM. 
Among the two most studied scaling relations are the radial profiles of temperature and total 
mass. The mass profile in particular is named the NFW profile after the authors of one of the 
original papers on this subject (Navarro et al. 1996). Deviations from the expected scaling in 
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high-quality data indicate the importance of non gravitational processes and/or unaccounted 
bias in the data. The temperature profile predicted by numerical simulations shows a significant 
drop with radius to about one third of the peak value at the virial radius (e.g. Loken et al. 
2002; Komatsu & Seljak 2002; Borgani et al. 2004; Roncarelli et al. 2006). Observationally, the 
temperature profiles are the key factor in deriving the cluster mass profile up to the virial radius. 
The precise mass profile will allow us to judge the validity of the present CDM framework, and 
gives assurance for our application of cluster properties to cosmological studies. 

For all of the above reasons it is important to extend X-ray cluster temperature mea- 
surements beyond the current limit of ~ 0.5 of the virial radius, particularly for relaxed clusters. 
Suzaku observed several such clusters in these regions, and some of the results have been pub- 
lished including PKS0745-191 (George et al. 2008), A1795 (Bautz et al. 2009), and A2204 
(Reiprich et al. 2009). All these clusters show a systematic trend of the temperature drop- 
ping to about one third of the central value, broadly consistent with theoretical expectations. 
However the statistical quality of the data for any individual cluster is limited and we need to 
look at many others to discern general behavior. The main difficulty for the measurement of 
ICM properties in the virial region is the low cluster surface brightness, which means that in 
no energy range does the cluster emission exceed the Galactic foreground plus cosmic X-ray 
background emission. Careful study of systematic errors is therefore mandatory when trying 
to assess the ICM properties around the virial radius. 

We have made Suzaku observations of A1413, a moderately distant cluster at redshift 
z = 0.1427 (Bohringer et al. 2000) whose size is well suited to our field of view. Assuming a 
Hubble constant of 70 km s~^ Mpc~^ or h-iQ = 1 as well as cosmological parameters of Q^q = 0.28 
and VIao = 0.72, we imply an angular diameter distance of 519 h^^ Mpc, a luminosity distance 
of 678 hjQ Mpc and a scale of 151.2 h^^ kpc per arcmin. 

Although we will measure this cluster's properties at the virial radius, we must make 
rather coarse spatial bins to do so. Thus we can not measure the actual virial radius with much 
precision. As a point of reference we adopt an often used nominal value, r2oo, which is the 
radius within which the cluster average density is 200 times the critical density needed to halt 
the expansion of the universe. For our cosmology 



in which (kT) is the cluster average temperature (Henry et al. 2009). An overdensity of 200 is 
contained within the virialized region of a spherical collapse in an Einstein-de Sitter universe 
at all red-shifts. Generahzing to a spherical collapse for our adopted cosmology at the red-shift 
of A1413, gives an overdensity of 109 for the viriahzed region (Henry 2000). However riog is 
only 22% larger than r2oo- So for comparison with previous work we adopt the latter as the 
nominal virial radius. 

The average temperature of A1413 integrated over the radial range of 70 kpc to rsoo is 




(1) 
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7.38±0.11 keV (Vikhlinin et al. 2006), where rsoo is defined analogously to r2oo, implying r2oo = 
2.24 hjQ Mpc or 14'. 8. Previous observations indicate the cluster is relaxed and there are high 
quality temperature and mass radial profiles available from both XMM-Newton and Chandra 
(Pointecouteau et al. 2005; Vikhlinin et al. 2006). There is some disagreement about the mass 
profile of A1413 in these two works. Pointecouteau et al. (2005) find rsoo = (1.13 ±0.03) /ifo^ 
Mpc and M500 = (4.82±0.42) x 10^^ h^f^ M©, while Vikhlinin et al. (2006) find (1.34±0.04) hjf^ 
Mpc and (7.79 ±0.78) x 10"'^'^ hjQ Mq, respectively, where M500 is the mass within rsoo- Note 
that both observations measure the temperature out to rsoo so the disagreement is not due to 
uncertainties in extrapolation. 

Throughout this paper, errors are at 90% confidence for one interesting parameter oth- 
erwise noted. 

2. Observations 

2.1. Suzaku 

We observed the northern region of A1413 with the Suzaku XIS detectors. In table 1, 
we give the details of our observation, and in figure 1(a), we show the XIS field of view (FOV) 
superimposed on the XMM-Newton image of A1413. The XIS instrument consists of 4 CCD 
chips; one back-illuminated (BI: XISl) and three front-illuminated (FI: XISO, XIS2, XIS3), with 
each is combined with an X-ray telescope (XRT). The IR/UV blocking filters had accumulated 
a significant contamination by the time of the observation since its launch (July 2005); we 
include its effects on the effective area in our analysis. The XIS was operated with normal 
clocking mode, in 5 x 5 or 3 x 3 editing modes. The spaced-row charge injection (SCI) was not 
applied, and all the four CCDs were working at the time of the observation. 

We show the FI±BI image in the 0.5-5 keV energy band in figure 1(b). The non X-ray 
background (NXB), cosmic X-ray background (CXB), and the Galactic background compo- 
nents (GAL) are subtracted as described below, and the result smoothed by a 2-dimensional 
gaussian with a — 16" are shown. The image is corrected for exposure time variations, but not 
for vignetting. Screening requirements are C0R2 > 8 GV and 100 < PINUD < 300 cts s~^ 
where C0R2 is the cut-off-rigidity calculated with the most recent geomagnetic coordinates 
and PINUD is the count rate from the upper level discriminatory of the Hard X-ray Detector 
(HXD) PIN sihcon diode detectors (see Tawa et al. 2008). The circles with 70" and 125" radii 
enclose excluded point sources. The small white circles indicate point sources detected in the 
XMM-Newton data. Blue circles show sources selected by eye in the Suzaku image. 

We used HEAsoft ver 6.4.1 and CALDB 2008-06-21 for all the Suzaku analysis presented 
here. We extracted pulse-height spectra in five annular regions from the XIS event files. The 
inner and outer radii of the regions were 2'. 7 — 7', 7' — 10', 10' — 15', 15' — 20', and 20' — 26', 
respectively, measured from the XMM-Newton surface brightness peak of A1413 at (R.A., Dec.) 
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Table 1. Log of SuzaJcu observations of Abell 1413 



Observation ID 

Date of observation 

Exposure time (ks) 

(C0R2 > GV) 

(C0R2 > 8 GV) 

(C0R2 > 8 GV and 100 < PINUD < 300 cts s-^) 

(R.A., Dec.) in J2000 * 

XIS mode 

^'n 



800001010 

2005-NOV-15 19:54:46 - 2005-Nov-18 14:14:45 

XIS0:107.4, XIS1:108.0, XIS2:107.5, XIS3:107.6 
XISO: 76.1, XISl: 76.4, XIS2: 76.2, XIS3: 76.2 
XISO: 71.9, XISl: 72.0, XIS2: 72.0, XIS3: 72.0 
(ll'^55'"19?0, 23°24'30") 

5x5/3x3, normal clocking, window off, SCI off 

2.19 X 1()2" ciir- (Dickry .V Lockiiiau 1!)!)()) 



* Average pointing direction of the XIS, shown by the RAJNTOM and DECNOM keywords of the FITS event files. 



= (ll'^SS^'lS.'^?, 23°01'48") in J2000. Wc analyzed the spectra in the 0.5-10 keV range for the 
FI detectors and 0.4-10 keV for the BI detector. In the the 2'. 7 — 7' annulus, we ignored the 
energy band 5-7 keV for the FI detectors when we analyzed the spectra, because those data 
were affected by Mn-K^^ (5.9 keV) X-rays from the ^^Fe calibration source. In other annuli, 
positions of the calibration sources themselves were masked out using the calmask calibration 
database (CALDB) file. 

2.2. XMM-Newton 

We analyzed an image in the energy band 0.35-1.25 keV taken with XMM-Newton 
(Pratt & Arnaud 2002). This observation was carried out in June 2000 (OBSID: 0112230501). 
The exposure time was 25.7 ks (M0S1,M0S2). SAS ver 6.0 and HEAsoft ver 6.4.1 were used 
for the analysis. XMM-Newton has much higher spatial resolution compared to Suzaku. We 
used this image as input for the response simulators and to find point sources. Pratt & Arnaud 
(2002) derived a ratio of minor to major axis to be 0.71 and a position angle 2°26' based on 
the XMM data. Since Suzaku coverage is limited in the north section of the cluster as shown 
in figure 1, we did not include the cluster ellipticity in our analysis. 

3. Background Analysis 

Accurate estimation of the background is particularly important when constraining the 
ICM surface brightness and temperature in the outer region of clusters. We assumed that the 
background is comprised of three components: non- X-ray background (NXB), cosmic X-ray 
background (CXB) and Galactic emission (GAL), which itself is comprised of two components. 
In this section we describe how we estimate all these background components. 



5 



(a) 



(b) 



+23°40' 



1 iV I •10- 

cts Ms"'' pixel"'' \j ,--•'* 



+23°30' : 



+23°20' 



+23° 10' 



■1"-"A1413 (0.35 - 1.25 keVf) 



+23°50' 



+23°40' 



+23°30' 



^^^B 

-2 -1 X 






10 ^"/V 




cfeMs'Vixer' f G* 


^ \ 
















/ / 

1 1 ^ 


" A1413{0 


^^^keV) ' 



Fig. 1. (a) XMM-Newton MOSl + M0S2 image of A1413 in tiie 0.35-1.25 keV band. The image is 
corrected for exposure, vignetting and background. The white and blue boxes show the fields of view of 
the Suzaku XIS and Chandra ACIS (Vikhlinin et al. 2006). The green circle shows r2oo of 14'. 8. Color 
scale unit is cts Ms~^ pixel"^ (1 pixel — 2". 49 x 2". 49). (b) Background subtracted Suzaku FI+BI image 
of the outskirts of A1413 in the 0.5-5 keV band smoothed by a 2-dimensional gaussian with a — 16". 
The image is corrected for exposure time but not for vignetting. Color scale unit is cts Ms^^ pixel^^ (1 
pixel = 1".04 X 1".04). C0R2 > 8 GV and 100 < PINUD < 300 cts s'^ screening was applied. The ^^Fe 
calibration source regions are also included in the figure, because they have negligible counts in this energy 
band. Large white circles denote 7', 10', 15', and 20' from the surface brightness peak of the XMM-Newton 
image. Small white and blue circles show the excluded point sources. 

3.1. Point Source Analysis 

We want to excise point sources because we are only interested in this paper in the 
ICM. However, since the CXB is comprised of faint point sources, we then need to correct the 
background level for the resolved sources. This and the next section describe the procedure we 
used for these tasks. 

We used the XMM-Newton image to detect point sources in the XIS FOV because its 
spatial resolution (14" half power diameter; HPD) is better than Suzaku's (2' HPD). We 
detected 10 point sources using wavdetect of CIAO, and extracted source and background 
spectra by setting the extraction radius of 33" and 33" — 66", respectively. First, we checked 
that the MOSl and M0S2 spectra of each source were consistent. Then we summed the 
MOSl and M0S2 spectra to increase the statistics, and fitted the spectrum of each source 
to evaluate individual spectral parameters. Finally we added the spectra of all the point 
sources to estimate how much of the CXB these sources resolve. We fitted the spectra by 
wabs X pegpwrlw. The best-fit parameters for the individual point sources and their sum are 
shown in table 2. We obtained x^/ dof = 87.2/77 for the power-law fit to the combined spectrum 
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(figure 2(a)), indicating a reasonable spectral fit. The photon index is F = 1.92 ±0.09 and the 
flux is 3.23toil X lO-^-'^ erg cm'^ s-\ 

We also searched for point sources located outside of the XMM-Newton field with Suzaku, 
finding an additional 5 sources by eye. These sources all show statistical significance higher than 
3.9(7 against the brightness fiuctuation around individual sources. We performed spectral fits 
to all the point sources with Suzaku according to the following procedure. The source photons 
came from a circle of 40" radius with encircling annular background region of 40" — 100" radii. 
We selected the source regions so they did not overlap each other. These source and background 
areas could be slightly different among the detectors and sources due to filtering by the calmask 
regions and the presence of hot pixels. We added the FI spectra from XISO, XIS2, and XIS3 
detectors, and summed the BACKSCAL keyword in the FITS header, which correspond to 
the area of extraction region, Agr or A]^^. Then, we carried out spectral fits for the FI and BI 
spectra simultaneously using the same spectral model as before, first for the individual sources 
and then the sum of all the point sources. 

We show the best-fit parameters for the individual point sources and their sum in table 2 
except for sources 04, 05, and 08, because they were faint so that we could not estimate their 
background reasonably. Obtained fluxes of the sources are shghtly affected by leaked photons 
of the target to the surrounding background regions. To correct for this effect, we calculated 
the ratio, /leak, of the leaked photons in each background region to the detected photons in 
the source region using the '^xissim" FTOOL (Ishisaki et al. 2007). We corrected the original 
source flux by multiplying a factor 1/(1 — /leak ^sr/^bg) — 1/(1 ~ 0-2 /leak) in the Fx columns of 
Suzaku in table 2. Figures 2(b) and (c) show the combined spectra of all sources for FI and BI. 
We obtained x^/dof = 113.1/117 for the power-law flt to the combined spectrum, indicating a 
reasonable spectral flt, too. The photon index is F = 1.82 ±0.12 and the flux is 4.83lo:56 x 10~^^ 
erg cm"^ s~^ (2-10 keV). 

The number of sources we found and their total flux are consistent with that expected 
from the logA^-logS* relation summarized in flgure 20 of Kushino et al. (2002). The detected 
sources ranges from ~ 10^^'^ to ~ 10^^^ erg cm~^ s~^. We excised all the point sources detected 
in either the Suzaku or XMM-Newton observations. Normally we excluded a region of 70" 
radius but used 125" radius for two sources (09 and 14 in table 2). 

3.2. Cosmic X-ray Background 

An ICM temperature measurement in the outer regions of a cluster is very sensitive to the 
CXB level. We took the 100% CXB surface brightness to be /q = 6.38 x 10"^ erg cm-^ s"^ sr-^ 
based on the ASCA-GIS measurements (Kushino et al. 2002). Moretti et al. (2008) summarized 
measurements (Gruber et al. 1999; McCammon et al. 1983; Gendreau et al. 1995; Vecchi et 
al. 1999; Kushino et al. 2002; Revnivtsev et al. 2003; DeLuca et al. 2004; Revnivtsev et al. 
2005; Hickox et al. 2006) of the CXB level, including their new result with XMM-Newton. The 
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Table 2. Best-fit parameters of detected point sources. 



Source 






FI+BI (S 


iizaku) * 




M0S1+M0S2 (XMM-Newton) t 


ID 


(R.A., Dec. 


) in J2000 


r 


F i 




/leak 


r 




xVdof 


01 


;il'^55'"38f7, 


23°34'02") 


^■^-0.8 


<2.4 


28.0/18 


1.53 


1 q+O-5 


^•^-1.2 


16.6/17 


02 


[11^55'°30?8, 


23°38'09") 


1.7 (fixed) 


<2.1 


14.7/12 


1.41 


1 q+0(' 

-■-•^-0.5 


9 y+2.3 


5.7/12 


03 


[;ilh55'°24?9, 


23°37'00") 


9 9+0.5 
^•^-0.4 


1 6+°-^ 


37.4/30 


1.34 


1 7+0.6 
' -0.5 




11.2/9 


04 


[1 1^55^21 ?6, 


23°33'02") 










-1 Q+0-3 




11.1/24 


05 


[llh55'°18?2, 


23°28'10") 










1 1+°-^ 


5.9l|9 


63.6/58 


06 


[llh55'°17?3, 


23°35'47") 


1 4+0.4 




28.1/27 


1.39 


^•^-0.2 




28.2/28 


07 


[llh55'°04?6, 


23°31'11") 


1 Q+0-3 
^■^-0.3 




39.2/37 


1.26 


9 1+0.2 
^•■'■-0.2 




40.9/37 


08 


[;ilh54™56?9, 


23°36'52") 










1 Q+0.5 
^■•^-0.5 


6.8ttt 


5.4/5 


09 


[;ilh54™51?7, 


23°34'49") 


1 7+0.2 
-0.1 


16.0t|^ 


79.1/58 


1.43 


9 9+0.2 
^•'^-0.2 




42.5/41 


10 


[llh54'°58?l, 


23°35'23") 


1.7 (fixed) 


< 1.8 


18.0/17 


1.21 


1 7+1.3 
-■^•'-1.2 


1 1+2-9 


5.1/7 


11 


^11^54'°33?0, 


23°37'51") 


1 7+0.6 
^■'-0.5 


q -1+2.3 


13.5/13 


1.24 








12 


[11^55™16?5, 


23°38'37") 






15.7/19 


1.44 








13 


[;il''55'"20?3, 


23°40'32") 


1+1° 


<2.8 


18.7/18 


1.44 








14 


[ll'>55'"29?7, 


23°47'26") 




2 7+1-'^ 

^•'-1.3 


22.4/14 


1.28 








15 


;ilh54'"47?7, 


23°39'23") 




6.6^?:^ 


23.2/19 


1.43 








Total 






1 00+0.12 
■'■•°^-0.12 


48.31^:° 


113.1/117 


1.36 


1 09+0.09 
-■-•^^-0.09 


32.3it^ 


87.2/77 



* Source-04, 05, and 08 are excluded because they exhibited negative counts after the background subtraction, 
t Source-11, 12, 13, 14, and 15 are out of MOSl and MOS2 FOVs. 
i Unit of flux is 10'^* erg cm-^ s"! (2-10 keV). 

measured CXB surface brightnesses show a significant range from the HEAOl value of (5.41 ± 
0.56) X 10-^ erg cm-^ g-^ sr-^ (Gruber et al. 1999) to (7.71 ± 0.33) x 10"^ erg cm-^ g-^ sr-^ 
with SAX-MECS (Vecchi et al. 1999) in the 2-10 keV band. Recent measurements show the 
flux to be within about 10% of the level reported by Kushino et al. (2002). 

We estimated the remaining CXB surface brightness after the above point-source sub- 
traction by the following three methods: (1) subtracting the summed point source fluxes mea- 
sured with Suzaku from the 100% CXB, (2) subtracting the summed point source fluxes esti- 
mated using the logA'^-logjS' relation, and (3) fitting a power-law model to the diffuse emission 
in the 20' — 26' region after the point sources arc excised. 

In case (1), we subtracted contribution of the excised sources of 1.80]'^o;2i x 10~* 
erg cm-2 s'^ sr'^ from the 100% CXB, dividing Fx = 4:.83tofQ x lO'^^ erg cm'^ s^^ of the 
Suzaku total by 17'. 8 x 17'. 8 area of the XIS FOV. In case (2), we calculated the integrated 
point source flux per steradian from 
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(a) M0S1+M0S2 



(b) FI (Suzaku) 



(c) BI (Suzaku) 
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Fig. 2. Power-law model fit to the sum of all point source spectra, (a) M0S1+M0S2, (b) FI, and (c) 
BI (black: source spectra, grey: best-fit model). 



I{S > So 



kn 



C-7+2 



7- 



(2) 



where k and 7 are the differential logA^-logS* normahzation and slope, respectively. We take 



1.58 X 10-15 sr-i 



[erg cm 
,-2 



-2 s-M^-i and 
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2.5, from Kushino et al. 



nominal values, ko 

(2002). 5*0 is taken as 2 x 10^^^ erg cm^^ s^^, which is slightly higher than our flux limit, 
because the assumed logA^-logS' in equation (2) does not take into account the flattening of 
the relation in the fainter flux end. In case (3), we flt the spectra from the solid angle in the 
20' — 26' annulus that remain after the source excision by a power-law model using a uniform 
flux ancillary response flle (ARF; see section 4.1). The ARF assumes that X-ray photons comes 
into the detectors uniformly from the sky direction within 20' radius from the optical axes of 
the respective XRTs. The model flt is apec+ wabs x [apec + powerlaw) where the two apec 
components represent the galactic emission. This is the 2T-III model described in section 3.4. 
In this case, value of the Jq — /x column is determined by the spectral flt, and then /x is 
calculated assuming Iq = 6.38 x 10^^ erg cm^^ st~^ in table 3. 

We summarize our estimations of the remaining CXB surface brightness, Iq — Ix, in 
table 3. All three methods give consistent results. Hereafter we will use a nominal diffuse 
cosmic X-ray background spectrum (after subtraction of point sources brighter than ~ 1 x 10"^^ 
erg cm-2 in 2-10 keV band) described by a power-law with a photon index F = 1.37, and 
surface brightness 4.73 x 10"^ erg cm-^ s"^ sr^^ in the 2-10 keV band, which comes from the 
2T-III (a) row of the method (3). We adopt this method because it directly measures the 
quantity of interest in our observations. 

To estimate the amplitude of the CXB fluctuations, we scaled the measured fluctuations 
from Ginga (Hayashida 1989) to our flux limit and FOV area. The fluctuation width is given 
by the following relation, 

0.5 / \ 0.25 

(3) 



'^Suzaku CTQjjjga^ / fio, Suzaku 



c, Suzaku 



-^CXB -^CXB \ ^e.Ginga / \ 'S'c,Ginga / 

where (o"suzaku/-^cxB) means the fractional CXB fluctuation width due to the statistical flue 



9 



Table 3. Estimation of the CXB surface brightness after the point source excision. 



(1) 4.58t°J2 1.80t°J2 1.41 (fixed) 

(2) § 4.15 2.23 1.41 (fixed) 

(3) 2T-III (a) II 4.73^0^1 LeslES;^^ 1-371^^ 

(3) 2T-III (b) » 4.69j:°;i« ^-(^^-oll l-40to:o7 

(3) 2T-III (c) II b.iet'oil 1.22l°i^ l.Ut'of, 

(3) contami+20% H 5-04tgi^ l-34tgi^ l-45lg:gi 

(3) contami-20% H 4.95^^]^ ^-^^-o fs l-44lg;g^ 



* Estimated surface brightness of the CXB after the point source excision in unit of 10~* erg cm~^ sr~^ (2-10 keV). 
t Contribution of the resolved point sources in unit of 10~* erg cm~^ sr~^ (2-10 keV). 
^ Assumed or estimated photon index of the CXB. 

§ Surface brightness of 100% of CXB is assumed as lo = 6.38 x 10~^ erg cm~^ s~^ sr~^ (2-10 keV). Integrated point source 

contribution, 7x, is calculated with equation (2). See section 3.2 for details. 

II See section 3.4 for definition. Abundance model is by Anders & Grevesse (1989). 

' See section 3.4 for definition. Abundance model is by Feldman (1992). 

tuation of discrete source number in the FOV. Here, we adopt CTGinga/^cxB = 5%, with Sc 
(Ginga: 6 x 10~^^ erg cm~^ s~^) representing the upper cut-off of the source flux, and 
(Ginga: 1.2 deg^) representing the effective beam size (or effective sohd angle) of the detector. 
We show the result, a/IcxB, for each spatial region in table 4. 

3.3. Non X-ray Background 

The non X-ray background (NXB) spectra were estimated from the Suzaku database of 
dark earth observations using the procedure of Tawa et al. (2008). We accumulated data for 
the same detector area, for the same distribution of C0R2 as the A1413 observation using the 
xisnxbgen FTOOLS covering 30 days before to 90 days after the observation period of A1413. 
To increase the A1413 signal-to-noise ratio by reducing the NXB count rate, we required C0R2 
to be > 8 GV and PINUD to be between 100 and 300 cts s~^. After this screening the exposure 
time dropped from 108 ks to 72 ks, nevertheless the fit residuals were reduced. We also tested 
other screening criteria, such as C0R2 > 8 GV and C0R2 > 5 GV, both with no PINUD 
screening. The former criterion did not affect the final spectral results significantly, but the 
latter gave different ICM temperatures. To test a possible NXB uncertainity systematic error, 
we varied its intensity by ±3% as investigated by Tawa et al. (2008). 

3.4. Galactic Components 

We fit the data in the 20' — 26' region to constrain the foreground Galactic emission, using 
the same uniform-sky ARF as the CXB component. We investigated the best model to use and 
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Table 4. Properties of the spatial regions used. 



Region * 


fie ^ 


Coverage 


SOURCE, t 


cr//cxB ^ 




FI counts (0.5-10 keV) " 






(arcmin 




RATIO. REG 


(%) 


OBS 


NXB 


CXB 


GAL 


/bgd 


2'.7-7' . 


. . 18.6 


14.2% 


2.60% 


15.4 


3, 828 ±62 


855 ±86 


560 ± 86 


81 ±9 


39.1 ±3.5% 


7' - 10' . 


. . 25.6 


16.0% 


1.49% 


13.1 


3, 568 ±60 


1,241 ±124 


966 ±127 


131 ±11 


65.5 ±5.3% 


10' - 15' 


. . 55.0 


14.0% 


1.40% 


9.0 


6, 340 ±80 


2, 428 ±242 


2, 460 ±220 


296 ± 17 


81.8 ±5.3% 


15' - 20' 


. . 86.5 


15.7% 


1.27% 


7.1 


9, 156 ±96 


4, 162 ±416 


4, 035 ±288 


499 ± 22 


95.0 ±5.6% 


20' - 26' 


. . 38.6 


4.5% 


0.47% 


10.7 


4, 547 ±67 


2, 523 ±252 


1,907 ±204 


272 ± 16 


103.4 ±7.3% 


Region * 


fie ^ 


Coverage 


t SOURCE. * 


cr//cxB ^ 




BI counts (0.4-10.0 keV) ^ 






(arcmin 




RATIO. REG 


(%) 


OBS 


NXB 


CXB 


GAL 


/bgd 


2'.7-7' . 


. . 18.4 


14.0% 


2.60% 


15.5 


2, 042 ±45 


748 ± 75 


208 ± 32 


85 ±9 


50.9 ±4.6% 


7' - 10' 


. . 25.5 


15.9% 


1.56% 


13.1 


2, 546 ±50 


1,088 ±109 


392 ± 51 


144 ± 12 


63.8 ±5.1% 


10' - 15' 


. . 54.9 


14.0% 


1.41% 


9.0 


4, 447 ±67 


2, 277 ±228 


1,012±91 


331 ± 18 


81.4 ±5.7% 


15' - 20' 


. . 87.4 


15.9% 


1.32% 


7.1 


7, 425 ±86 


4,418±441 


1,857 ±132 


631 ±25 


93.0 ±6.3% 


20' - 26' 


. . 24.6 


2.8% 


0.37% 


13.3 


2, 984 ±55 


1,954 ±195 


706 ± 94 


264 ±16 


98.0 ±7.5% 



* Radii are from the XMM-Newton surface brightness peak in figure 1(a). 
t The average value of the four detectors. 

t SOURCE.RATIO.REG= Coverage x JJ^"" S(r)rdr/ S (r) r dr , where S{r) represents the assumed radial profile of A1413. 

We confined S{r) to a 60' x 60' region on the sky. 

§ 5c = 1 X 10~^* erg cm~^ is assumed for all regions. 

II OBS denotes the total observed counts. NXB, CXB and GAL are the estimated counts, /bgd = (NXB + CXB + GAL)/OBS. 

the best-fit model parameters. In all cases, we also included a power-law model to represent 
the CXB. We tried a single temperature thermal plasma model, IT: apec + wabs x powerlaw, 
a two temperature model, 2T: wabs x [apec^ + apec2 + powerlaw), and a two temperature 
model following Tawa et al. (2009), 2T-III: apeci + wabs x (^apec2 -\- powerlaw) . In all models, 
redshift and abundance of the apec components were fixed at 0.0 and 1.0, respectively. The 
two temperature variants try to model the Local Hot Babble (LHB) and the Milky Way Halo 
(MWH). We tried three types of the 2T model: both temperatures fixed to 0.204 keV and 
0.074 keV given by Lumb et al. (2002), one temperature fixed to 0.074 keV and the second 
temperature free, both temperatures free. We call the first model as 2T-I, and the second model 
as 2T-II. The third model did not converge in the fitting, so that we do not discuss it further. 

We found that the IT and 2T-I models gave worse values compared with the 2T-II 
and 2T-III fits. We show the best-fit parameters in table 5 for the 2T-III model, which we 
adopt. We find that the LHB and MWH temperatures are 0.112tg:gg| keV and 0.278^^;^?^ keV, 
respectively. These values are consistent with those obtained by Tawa et al. (2009). We also 
show in table 5 how the best-fit parameters change as a result of systematic changes in the 
CXB and NXB levels and of the abundance model used (labeled (a) or (b)). The variations 
are small: less than ±10% for the temperatures and ±15% for the normahzations. Finally, 
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Table 5. Galactic components best fit parameters and 90% confidence errors. 



kTi (keV) Normi § Si H /cTa (keV) Norm,2 § ^2 " 

2T-III (a) * 0.112t°:™^ 1.994^°;!^^ O.eie^^^^; 0-^78troll OAMI^oIt 0.2581^:! 

2T-III (b) t 0.1101°:^°^ 2.199tVitl Om^trotl 0.314lO:S^« 0.222iS;°l^ 0.3021°:! 

2T-III (c) ^ 0.113t°;°°^ 1.727tri^l OMOt^oll 0.260^^0^ 0.20i;°;°l 0.2661°;: 

NXB+3%,CXBmax ... 0.112t°:EJ°-^ 2.015l°:l«^ 0.640tS:°^^ 0.31ltJ;;,g 0.197tg;g^^ 0.2691°;! 

NXB-3%,CXBmin .... 0.11llO:2;,^6 2.170l°J3^ 0.6511^^? 0.319i°;™^ 0.227lS;^?3 Qgi^+o. 

contami+20% 0.11ll°°;2 2.254l°i27 0.660tg:g^^ 0.269lg;^l| 0.223tg:g|^ 0.295lg;l 

contami-20% 0.113tH^^ 1.79lt°n]f^ 0.569t°afA 0.286tn?a 0.173lH'4 0.232t0-: 



* Abundance model is by Anders & Grevesse (1989). 
t Abundance model is by Feldman (1992). 

t Including two gaussian models of Ovii and Oviii- Abundance model is by Anders & Grevesse (1989). 
§ Normalization of the apec component scaled with a factor of SOURCE^RATIO^REG /^e in table 4, 

Norm = SOURCE.RATIO.REG/ne J neUndV / {■in {1 + z)'^ D^) x IQ-^o cm-^ arcmin-^, where Da is the angular diameter 
distance to the source. 

II Surface brightness in unit of 10~® photons cm~-^ arcmin"^ (0.5-10 keV). 



our baseline CXB+GAL model is denoted 2T-III (a), apeci + wabs * {apec2 + power law) with 
abundances from Anders & Grevesse (1989). We link all parameters of this model, except an 
overall normalization, when performing the fits for the different spatial regions described in 
section 4. 

3.5. Background Fraction in Each Region 

Table 4 presents many properties of the spatial regions we analyzed. The columns are 
the annular boundaries; the actual solid angle of each region observed, Qe! the coverage fraction 
of each annulus which is the ratio of Jig to the total solid angle of the annulus. Coverage; the 
fraction of the simulated cluster photons which fall in the region compared with the total 
photons from the entire simulated cluster, SOURCE_RATIO-REG; the CXB fluctuations due 
to unresolved point sources, c/Icxb', the observed counts, OBS; the estimated counts for each 
background component, NXB, CXB, and GAL; and the fraction of background photons given 
by /bgd = (NXB+CXB+GAL)/OBS. 

The NXB count rates are calculated from the dark earth data. We simulated the CXB 
and GAL components spectra using xissim with the flux and spectral parameters given in row 
2T-III (a) of tables 3 and 5, assuming a uniform surface brightness that fllls the fleld. We plot 
the NXB and CXB spectra in figures 3 and 4. These spectra gave the count rates in table 4. 
In the outermost region of 20' — 26', /bgd is consistent with 100%. This confirms the accuracy 
of our background estimation. 
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4. Spectral Analysis 



4-1- Spatial and spectral responses 

We need to prepare the spatial and spectral responses which are necessary for reducing 
and analyzing our observations of A1413. These responses have complicated properties for 
extended sources. Indeed they depend on the surface brightness distribution of the source and 
so are unique for each annular region. Monte Carlo simulators are used to generate some of 
the responses. The X-ray telescope + XIS simulator is called xissim, and the ARF generator 
using the simulator is called xissimarfgen (Ishisaki et al. 2007). We used version 2008-04-05 of 
the simulator. 

A surface brightness distribution is necessary for xissim and xissimarfgen, because the 
point spread function (PSF) of the XRT produces an efficiency that is correlated among adjacent 
spatial cells. Since the XIS FOV did not include the brightness peak of A1413, we used the KBB 
model of Pratt & Arnaud (2002) to generate the ARF. We numerically projected the KBB 
3-dimensional model of the gas density to generate the input surface brightness distribution. 
Since the ARF describes the detection efficiency as a function of energy, no particular spectral 
shape is required for input. The effect of the XIS IR/UV blocking filter contamination is 
included in the ARF based on the calibration of November 2006. The normalization of the 
ARF is such that the measured flux in a spectral fit for a given spatial region is the flux from 
the entire input surface brightness. The flux just from the spatial region is the fit fiux times 
the xissimarfgen output parameter SOURCE_RATIO-REG (table 4). The surface brightness 
from a given spatial region is then the usual flux from the region divided by the solid angle 
that subtends from the observer. 

We examined how many photons accumulated in the five annular regions actually came 
from somewhere else on the sky because of the extended telescope PSF. We show in table 6 the 
results for the FI+BI detectors in the 0.5-5 keV band. These numbers agree well within 1% 
for individual sensors and other reasonable energy bands. About 70% of the photons detected 
in each region actually come from the corresponding sky region. Serlemitsos et al. (2007) gives 
an upper limit on the error in the simulation at 20'. He reported that the actual stray intensity 
levels were less than twice those predicted by xissim due to the XRT reflector alignment errors 
and reflections from the pre-coUimator blades. 

The redistribution matrix flle (RMF), which gives the spectral response to a mono- 
energetic input, is the same for all sources. It was generated with xisrmfgen version 2007-05-14. 
Degradation of the energy resolution is included based on the calibration in November 2006. 

4.2. Model Fitting 

We used XSPEC version 12.4.0y for all spectral fitting. The FI and BI spectra were 
fitted simultaneously. We employed a wahs x apec model for the ICM emission of the cluster. 
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(a) FI, 2'.7-7' 



(b) FI, r - 10' 




Energy (keV) Energy (keV) 

(c) FI, 10' - 15' (d) FI, 15' - 20' 




0.5 1 2 5 10 - 0.5 1 2 5 10 



Energy (keV) Energy (keV) 

(e) FI, 20' - 26' 



O : 




0.5 1 2 5 10 



Energy (keV) 

Fig. 3. Spectra for the individual annular regions observed with the FI sensors. The total observed 
spectrum minus the estimated NXB is the black crosses, the estimated NXB is the grey crosses, and the 
fitted CXB component is the solid line. The screening used are C0R2 > 8 GV and 100 < PINUD < 
300 cts ^^Fe calibration source regions, namely calmask, are excluded except for (a). 
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(a) BI, 2'.7-7' 



(b) BI, r - 10' 
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Fig. 4. Same as figure 3, but for the BI detector. All the ^^Fe calibration source regions are excluded. 
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Table 6. Emission weighted radius and estimated fractions of the ICM photons accumulated in detector regions coming 
from each sky region for FI+BI in the 0.5-5 keV band. 



Detector Emission weighted 

region radius * 

2'.7-7' 4.7^2:^ 

7' - 10' 8.0ti:g 

10' -15' 

15' - 20' 18.6^3:6 

20' - 26' 



Sky region 

0-2.7' 2.7-7' 7-10' 10-15' 15-20' 20-26' 

21.5% 73.2% 5.1% 0.2% 0.0% 0.0% 

16.0% 21.8% 54.9% 7.3% 0.1% 0.0% 

6.7% 7.3% 14.0% 67.3% 4.7% 0.0% 

4.1% 2.7% 2.7% 16.8% 67.1% 5.4% 

6.7% 6.7% 0.0% 0.0% 20.0% 66.7% 



* Emission weighted radius from the surface brightness peak of the XMM-Newton. 

The wabs component models the photoelectric absorption by the Milky Way, parameterized 
by the hydrogen column density that was fixed at the 21 cm value (Dickey & Lockman 1990). 
The apec is a thermal plasma model. Its fitting parameters are normalization, kT and the ICM 
abundance. The redshift was fixed at the optical spectroscopic value {z = 0.1427). Additional 
fitting parameters are the two normalizations and temperatures of the GAL components, and 
the normalization and photon index of the power-law model for the CXB component, as de- 
scribed previously. We did not fit the ICM component in the outermost 20' — 26' region because 
we can explain the observed spectrum without it, as we show in figure 5(e). This situation was 
planned as we wanted to have an in-field measurement of the background. In figures 3 and 4, 
we compare the intensities of the observed spectra minus the NXB to the spectra of the NXB 
and CXB components. Figure 3(a) shows very strong Mn-Ko, fine at 5.9 keV from the ^^Fe 
calibration source, therefore we ignored the 5-7 keV energy band when we fit the FI spectrum 
of this annulus. 

4-3. Results 

In figure 5, we show the best-fit spectra in each spatial region. These figures show the 
observed spectra after subtraction of the NXB, as well as the best-fit. These figures show 
that individual spectra are well fitted by the model in each region. The normalization for 
the ICM component was fixed to zero in the 20' — 26' annulus to estimate the background. 
The ICM spectra did not show strong emission lines. Because of the low S/N ratio, it was 
difficult to constrain the model parameters in the 15' — 20' annulus. Therefore, we finked the 
ICM temperature and abundance in this region to that of the region next interior to it, the 
10' — 15' annulus. The best-fit parameters were consistent within the systematic errors for the 
two regions. The emission weighted average radius for the combined region is 12'.42^^/q7. 

Table 7 shows the best-fit parameters for the ICM model in each region. We fitted with 
two different solar abundances, namely Anders & Grevesse (1989) and Feldman (1992). The 
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Table 7. Best fitting parameters of the spectral fits with 90% confidence errors for one parameter. 



2T-III (a) * kT Abundance Norm § S " xVdof 

(keV) (Zq) 

2'.7-7' 7.031^;^^ 0-44to;39 ^^-^^-1^26 ^■'^'^-oaI 77.4/107 

7' -10' 4.13ir6^ 0.54toii 4.53l°i° 2.121°^^ 98.7/116 

10' -15' ^-QOtoil 0.39+°oll 2.29tg:^^ 0.90to?o 130.1/118 

15' -20' t t 0.82tg:^^ OMtoAo 109.5/116 

20' -26' — — — — 152.7/113 

Total — — — — 568.4/570 

2T-III (b) t kT Abundance Norm § S H xVdof 

(keV) (Zq) 

2'.7-7' 7.Utlf^ 0.58lHo 16-04t^:l7 5.75t°:|^ 77.1/107 

7' -10' 4.41^0;?^ 0.66t° j^ 4.43^°^^ 2.11^°;^^ 100.6/116 

10' -15' 4:mtoil 0-77to:5i 2.07tg:J^ 0.90tg;g^ 129.6/118 

15' -20' t t 0.72tg;^^ 0.3ltg;5o 114.7/116 

20' -26' — — — — 149.2/113 

Total .",71.3/570 

2T-III (c) i kT Abundance Norm § S H xVdof 

(keV) {Zq) 

2.7' -7' 7.20ti:^^ 0.431^1? 26.54^^:^^ 9.4ltg;^^ 76.7/105 

7' -10' 4.33tg:?o 0-68to:22 ll-02to:60 5.44tg:^^ 99.8/114 

10' -15' 3.97tg:^6 0-53to:2i 2.07^^;" 0.89^^;!^^ 125.2/116 

15' -20' t t 0.66t^:}3 0-34to:o7 104.3/114 

20' -26' — — — — 154.5/113 

Total — — — — 560.5/562 



* Abundance model is Anders & Grevesse (1989). 
t Abundance model is Feldman (1992). 

^ Including two gaussian models of Oyii and Oyiii WHIM emission. Abundance model is Anders & Grevesse (1989) 
§ Normalization of the apec component scaled with a factor of SOURCE.RATIO-REG/Clg, in table 4, 

Norm = SOURCE.RATIO.REG/ne J UenndV / (47r (1 + z)'^D^) x IQ-^o cm-^ arcmin-^, where is the angular diameter 
distance to the source. 

II Surface brightness in unit of 10~® photons cm~^ arcmin"^ (0.4-10 keV). 



17 



Table 8. Same as table 7 except NXB±3%, CXBmax and CXBmin and contami±20%. Abundance model is Anders & 
Grevesse (1989). 



NXB-3%, CXBmin kT Abundance Norm * S t xVdof 

(keV) {Zq) 

2'.7-7' 7.57l:lil O-^Ttoif 16-94t?:^g 5.99tg;^| 78.7/107 

7' -10' 4:.84to,l\ O.eotg jg 4:Mtoil 2.34tg;^^ 98.6/116 

10' -15' 4.64t°:?? O.Sltgj^ l-07tg:l^ 0.43^^7 130.6/116 

15' -20' t t 0.98tg:l^ 0.4ltg;g^ 114.2/116 

20' -26' — — — — 157.1/115 

Total .",79.1/572 

NXB+3%, CXBmax kT Abundance Norm * S t xVdof 

(keV) (Zq) 

2'.7-7' Q-QOtlil 0-40lo.40 15.96t;:|| 5.5ltg;|^ 76.6/107 

7' -10' 3.59tg:^^ 0.53tg:25 4.16tg:^^ l-87to:32 104.3/116 

10' -15' 2.52t°j^ 0.35t°o\l 2.14l°;i^ 0.76^°;?? 130.3/116 

15' -20' t t 0.53tg:i^ 0.18tg;S^ 118.6/116 

20' -26' — — — — 150.1/115 

Total — — — — 579.9/572 

contami+20% kT Abundance Norm * S t x^/dof 

(keV) (Zq) 

2'.7-7' 6.89t}:^t 0.45+5];^;] 16.3ll}i^ 5.741°;^ 77.7/107 

7' -10' 4.011°:^^ 0.54t°;2^ 4.54l°:^2 2.10^°;^^ 99.0/116 

10' -15' 3.17tgji Q-'^Q-oaI 2.4lt°:g 0.901^0 131.3/118 

15' -20' t t 0.841^:^^ 0.30l^;2^ 109.6/116 

20' -26' — — — — 153.4/113 

Total — — — — 571.0/570 

contami-20% kT Abundance Norm * S ^ X^/dof 

(keV) {Zq) 

2'.7-7' 7.081};^^ 0.42+° 16.361?;^^ 5-78lo:|8 77.4/107 

7' -10' 4.19t°:^^ 0-54tl^;25 4.49l°:3i 2.nt°o fi 99.0/116 

10' -15' 3.82l°:^^ 0.44t°;^^ 2.211°;^} 0.89lJ]:2^ 128.4/118 

15' -20' t t 0-79;o:23 0.3ll^:^g 109.1/116 

20' -26' — — — — 153.2/113 

Total _ _ _ _ 567.1/570 



* Normalization of the apec component scaled with a factor of SOURCE JIATIO-REG /Q.c in tabic 4, 

Norm = SOURCE. RATIO. REG /Vie J neUndV / (47r (1 + z)^Dl) x IQ-^o cm-5 arcmin-2, where Da is the angular diameter 
distance to the source. 

t Surface brightness in unit of 10~® photons cm~^ s~^ arcmin"^ (0.4-10 keV). 



(a) FI+BI, 2'.7-7' 



(b) FI+BI, r - 10' 



2 10-3 



10-; 



< 



- 




1 






'l|lilltllllFf|l 


m 







0.5 

(c) FI+BI, 10'- 15' 

0.1 



10 



Energy (keV) 



> 

0} 



0.01 



2 10-3 



10-" 




0.5 1 2 

Energy (keV) 



10 



(e) FI+BI, 20' - 26' 

0.1 





Energy (keV) 



(d) FI+BI, 15' -20' 



0.1 



> 

CD 



0.01 



2 10-3 



10-" 



1 






: I.I 1 1 1 M 




(Kill 




1 ■ ■ -L 


1 


..|'Ti'Tijt' 

1 





0.5 



10 



Energy (keV) 



Energy (keV) 

Fig. 5. The upper panels show the observed spectra after subtracting the NXB, that is fitted with the 
ICM: wahs x apec model plus the GAL+CXB: apec^ + wahs x {apec2+ powerlaw) model in the energy range 
0.5-10 keV for FI and 0.4-10 keV for BI. The annular regions are: (a) 2'. 7 - 7', (b) 5' - 10', (c) 10' - 15', 
(d) 15' - 20', and (e) 20' - 26'. The symbols denote BI data (red crosses), FI data (black crosses), CXB 
of BI (purple), apeci of BI (grey), wahs x apecj of BI (light blue), ICM of BI (orange), the total model 
spectra of BI (green), and that of FI (blue). The lower panels show the residuals in units of a. 
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Fig. 6. Radial profiles for (a) temperature, (b) surface brightness (0.4-10 keV), (c) abundance, and (d) 
3-dimensional electron density. Red diamonds show our Suzaku results assuming the metal abundances of 
Anders & Grevesse (1989). Orange line indicates the best-fit profile using the Feldman (1992) abundances. 
Chandra results by Vikhlinin et al. (2005) are the black crosses, and the cyan crosses are the XMM-Newton 
results by Snowden et al. (2008). The uncertainty range due to the combined ±3% variation of the NXB 
level and the maximum/minimum fluctuation of CXB is shown by two blue dashed lines. We show by 
magenta dashed lines the uncertainties induced by a ±20% uncertainty in the amount of contamination 
in the IR/UV blocking filters. We also show in panel (b) the CXB level (horizontal dashed line) and the 
Galactic emission (horizontal solid line) . 
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derived abundance values are higher when we adopt the Feldman (1992) abundance, than the 
Anders & Grevesse (1989) case, because the Fe abundance relative to H in the former model is 
lower than the latter. 

In figure 6(a), we show temperature profiles observed with Chandra (Vikhlinin et al. 
2005), XMM-Newton (Snowden et al. 2008), and Suzaku (this work). These profiles are con- 
sistent with each other in the range 7' — 15'. The Chandra temperatures are about 20% higher 
than the XMM-Newton values at 2'. 7 — 7'. The tendency that Chandra gives higher tempera- 
ture than XMM-Newton typically becoming significant above kT ~ 5-6 keV is pointed out in 
figure 12 of Snowden et al. (2008). This discrepancy is due mainly to a Chandra calibration 
problem, namely the ground cahbration of the HRMA effective area had some errors espe- 
cially at the Ir edge (0.62 keV), and there also was uncertainity about the IR/UV blocking 
filter contamination. These uncertainties caused a large discrepancy between the Chandra and 
XMM-Newton measurements for high-temperature clusters. Recent updates of the Chandra 
CALDB, HRMA AXEFFA version NOOOS,^ corrected most of this discrepancy. However, there 
still remains some differences in cluster temperature by about 10% especially in hot objects. 
For temperatures below ~ 5 keV, Chandra and XMM-Newton results are mostly consistent 
with each other. 

We therefore used the XMM-Newton temperatures measured by Snowden et al. (2008). 
In fact, their values are higher than those of Pratt & Arnaud (2002) who used the same data set. 
This difference may partly be due to the different backgrounds used. Therefore, we assigned 
rather large errors of 10% even in the inner region of r < 2'. 7 for these data. We will quantify 
the systematic error of the Suzaku ICM temperature in the following section. 

We plot the related quantities, surface brightness, Sx, and 3-dimensional electron density, 
TT-e, in figures 6(b) and (d). We derived the Chandra surface brightness from the emission mea- 
sures provided by A. Vikhlinin (private communication). The XMM-Newton surface brightness 
is from Snowden et al. (2008). The Suzaku surface brightness comes from the normalization 
of the apec model fit. The surface brightness results are consistent with each other within 10'. 
In the outer region, the Suzaku surface brightness is significantly higher than the Chandra 
values. The cause of this discrepancy could be the different region of the cluster observed. In 
particular, Suzaku observed mainly along the major axis, while Chandra observed the minor 
axis, as we show in figure 1(a). We obtained the electron density by deprojecting the emission 
measure with method describe in Kriss et al. (1983). 

We show the abundance profile in figure 6(c). Our nominal values are higher than the 
results of Chandra and XMM-Newton. However, our errors are large and it is difficult to draw 
firm conclusions. 

^ http:// cxc.harvard.edu / ciao4. 1 / why/caldb4. 1 . 1 Jirma.html 
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4.4- Systematic Errors 

To estimate the systematic errors on our electron density, temperature and abundance 
profiles, we examined the effects of varying the background spectra from their nominal levels. 
We adopted a systematic error for the NXB intensity of ±3% and the level of the CXB fluctu- 
ation was scaled from the Ginga result (Hayashida 1989) as shown in table 4. We considered 
a ±20% error for the contamination thickness on the IR/UV blocking filters in front of the 
XIS sensors. As mentioned earlier, we also looked into the effect of the difference between the 
Anders & Grevesse (1989) and Feldman (1992) abundance models. 

We give the outcome of these variations in figure 6 and table 7 for the abundance model 
comparison, and in figure 6 and table 8 for the other comparisons. Systematic variations of the 
surface brightness are comparable to its statistical error for all the systematics we examined. 
The same is true of the temperature except for uncertainties on the UV/IR filter contamination, 
where the maximum possible range allowed is about 40% larger than the nominal statistical 
errors. Systematics on the abundance profile were less than the statistical uncertainties except 
for the outer two spatial bins with the Feldman (1992) abundance models. We conclude from 
this investigation that our statistical errors also encompass most possible systematic effects. 

4-5. Search for WHIM lines 

We searched for the warm-hot intergalactic medium (WHIM) which could exist in the 
filaments of large-scale structures of the universe. The outer regions of clusters may be con- 
nected to these filaments and are considered to be promising regions to search for possible 
WHIM emission. 

We analyzed the regions 2.7 - 7', 7' - 10', 10' - 15', and 15' - 20'. We fitted the FI+BI 
spectra simultaneously. We added two gaussian lines to model the oxygen emission lines. They 
had fixed redshifted energies of 0.508 keV (Oyii) and 0.569 keV (Oviii), with a fixed width of 
a = 0.0. The ICM spectra fitted with the additional two gaussian lines are shown in figure 7, 
and table 7(c) gives the fit results. The best temperatures are consistent with the results of 
the previous fit without the lines. Because redshifted line energies overlapped with those of the 
Galactic lines, we were unable to distinguish these emission lines directly. Table 9 gives our 
result for the line intensities which are either 2a upper limits or marginal detections. 

5. Discussion 

5.1. Temperature and brightness profiles 

Numerical simulations indicate that the intracluster gas is almost in hydrostatic equi- 
librium within the virial radius. For example, Roncarelli et al. (2006) showed that the radial 
density profiles are smooth out to ~ 2r200) while the electron temperature profile has a dis- 
continuity around 1.3-1.5 r2oo- Eke et al. (1998) performed hydrodynamic simulations in a 



22 



(a) FI+BI, 10' -15' 



(b) FI+BI, 15' -20' 




Energy (keV) Energy (keV) 

Fig. 7. Ovii (cyan) and Oyiii (pink) line spectra in 10' — 15' and 15' — 20' annuli. 

Table 9. Intensity of redshifted Oyil (0.508 keV) and Oyiii (0.569 keV) lines in unit of 10~^ photons cm~^ arcmin"^ 
with la upper hmits or 90% confidence errors for a single parameter. 



Region S'oviii ■S'ovii 

2'.7-7' < 0.119 < 0.135 

7' -10' < 0.075 < 0.091 

10' -15' < 0.085 0.094j:°:°^ 

15' -20' < 0.095 O.OSlt^ EJ^ 



ACDM universe, and discussed the possibility of nonequilibrium around rioo because the ratio 
of kinetic to thermal energy gradually increased from the center to this radius. 

Recent X-ray studies of the outer regions of clusters of galaxies with Chandra and XMM- 
Newton showed significant negative temperature gradients out to a typical radius of rsoo which 
is about half of r^m (Vikhlinin et al. 2006; Pratt & Arnaud 2002; Snowden et al. 2008). Even 
though the errors are large, it is significant that our temperatures continue this steady decline, 
going from about 7.5 keV near the center to ~ 3.5 keV at r2oo- Recent Suzaku results for the 
A2204 (Reiprich et al. 2009), PKS0745-191 (George et al. 2008), and A1795 (Bautz et al. 
2009) clusters also show a temperature drop to 2-3 keV at r2oo- The similar temperatures at 
r2oo are at least partly due to the fact that all these clusters have similar average temperatures 
of 5-7 keV. What is likely more significant is the factor of ~ 2 decrease in all cases. 

We attempted to compare our measured temperature and surface brightness profiles with 
theoretical predictions for relaxed clusters. Suto et al. (1998) gave ICM properties for clusters 
whose potentials follow NEW (Navarro et al. 1996) and modified NEW models, assuming that 
the ICM can be described by a polytrope. These models have 6 parameters and give a wide 
range of temperature and density distributions with radius. 
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We found that, although we could fit either one of the temperature or surface bright- 
ness profile with the model, it was not possible to fit both profiles simultaneously despite an 
exhaustive search of the 6-parameter space. When we fixed the scale radius to be = 350 kpc 
and jointly fit the temperature and brightness profiles, we obtained reduced values of 2.0 
using only the Chandra data and 3.7 for combined Chandra and Suzaku data, respectively. 
The likely reason for this result is that the ICM is out of equilibrium in the outer regions of 
the cluster. We examine this hypothesis in the next section using the entropy profile. 

5.2. Entropy profile 

Entropy carries information about the thermal history of the ICM, which is thought to 
be heated by accretion shocks outside the virial radius. The central regions of clusters often 
exhibit complicated physical phenomena, such as AGN heating and cooling flows, therefore it 
is difficult to trace the long-term evolution of clusters there. In contrast, the outer regions of 
clusters is where signatures of the structure formation history can be more clearly seen with 
the entropy profiles. We use the customary X-ray astronomy definition of entropy as 

S = kTnf/\ (4) 

We show the entropy profile derived from our data in figure 8(a). To compare the observed 
proflle with simulation results, we flt the XMM-Newton data from 0'.5 to 7' and the Suzaku 
data from 7' to 20' with a power-law model, given by 5" oc r'^. The XMM-Newton data outside 
of 7' have poorer quality than the Suzaku data, and one Suzaku point inside of 7' was also 
excluded because it is near the fleld edge with rather low data quality. 

We found the best-flt power-law indices to be 0.90 ±0.10 in 2' to T and 0.97 ±0.48 in 7' 
to 20'. The dividing radius of 7' corresponds to 0.47 r2oo- If we fit all the 7 data points from 2' 
to 20', then the slope becomes 0.90 ±0.12. These results indicate that there is no difference in 
the entropy slopes between the inner and outer regions. 

Voit (2005) reported S oc r^-^ based on numerical simulations of adiabatic cool gas ac- 
cretion, and our observational result shows a significantly flatter slope, at least for r < 7'. This 
feature is similar but less pronounced to those reported for A1795 (Bautz et al. 2009) in which 
the power-law index flattened (7 ^ 0.74) for r > 4' ~ 0.15 r2oo and for PKS0745— 191 where 
George et al. (2008) also found a flatter entropy profile in the outer regions. Our result for 
A1413 suggests that the entropy profile starts to fiatten from ~ 0.2 r2oo- To compare the en- 
tropy profiles with the simulated slope of 1.1, we divided the entropy by 5" oc r^'^ as shown in 
figure 8(b). There appears to be a deviation from the numerical simulation in the range of 
r > 0.2 r2oo, indicating the flattening of the entropy proflles. We note that the flattening is 
common to three clusters 

We compare our result with a hydrodynamical simulation by Takizawa (1998) that al- 
lowed for different electron and ion temperatures. We fit a /3-model density profile (parameters 
no,rc,/3) and a polytrope electron temperature profile (parameter polytrope index 7p) using the 
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simulated data in his tables 1 and 2. The resulting entropy profile shows a slope of 7p = 0.42 in 
the outer regions for the case of flat universe with (QqiAq) = (0.2,0.8). Even though this result 
might be an extreme case, it shows that a difference in the electron and ion temperatures can 
cause a flattening of the entropy profile. 

5.3. Equilibration timescale 

Ions carry most of the kinetic energy in the cluster outskirts, and they will be thermalized 
fairly quickly after accretion shocks or mergers. However, heating the electrons takes a long 
time because of the inefficient energy transfer between ions and electrons; the equilibration 
time for electron-ion coUisions (^ei) is about 2000 times longer than electron-electron process 
(tee) and about 45 times longer than ion- ion relaxation time (tn). 

According to Fox & Loeb (1997), Takizawa (1998), and Rudd & Nagai (2009), the 
electron-ion timescale including contributions from both protons and He^"^ is estimated as 
(Spitzer 1956) 

where In A is the Coulomb logarithm. We simply assume that ions are initially heated through 
accretion shocks at r2oo- In the post-shock region, ions achieve thermal equilibrium with a 
timescale of ta after this heating. The ion temperature Ti will then be significantly hotter than 
the electron temperature Tg. Eventually, thermal energy is transferred from ions to electrons 
through Coulomb collisions, and will equal Tj after the relaxation time tgi- 

We can compare the position-dependent time since the shock heating, tciapscd, with the 
equilibration timescale tgi- If ^ci is longer than tciapscdi then Tg would be expected to be signifi- 
cantly lower than Ti at that position. Denoting the velocity of inward propagation of the shock 
front as Wshock, we obtain 

'^200 — r ~ ^elapsed 1'shock- (6) 



The free-fall velocity of the gas at r2oo is ■i^fr,2oo = ^2GM2oo/ ^200- Using the strong shock 
approximation and neglecting the post-shock gas velocity compared with I'shock, Takizawa (1998) 
found 

t'shock - ^ t'ff,200- (7) 

Then, we can derive 

^ o ^200 — ^ /„s 
^elapsed — "J ) \°) 

which is independent of M2oo- In figure 8(c), we show teiapsed and tei- In the region outside of 
r ~ 0.9 r2oo, ^ei is significantly longer than teiapsed- Based on this calculation, it is likely that 
and Ti are significantly different in the outskirts of the A1413 cluster. 
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(a) Entropy (b) Entropy / r^^^ 



o 




Radius (r/r^J Radius {r/r^J 

Fig. 8. (a) Entropy profiles (black diamond: Suzaku, grey diamond: XMM-Newton, black solid line: 
fitted model to Suzaku in 7' — 20', black dashed line: fitted model to XMM-Newton in 0'.5 — 7', grey solid 

cross: PKS0745— 191, grey dotted cross: A1795). (b) Entropy normalized to ocr^-^ profile, (c) td profile 
(diamonds) compared with tdapsed (black solid line), (d) Te/Tgas profiles compared with the simulated 
result by Rudd & Nagai (2009). 
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5.4- Difference between Electron and Ion Temperatures 

Fox & Loeb (1997) were the first to investigate the two-temperature nature of the ICM. 
Takizawa (1998) showed that in a one dimension numerical simulation there existed a significant 
difference between the electron and ion temperatures, which will affect the entropy profile 
and the inferred gravitational mass. Recently, Rudd & Nagai (2009) reported the results of 
simulations which indicated that the temperature difference had a maximum of about 30% at 
'1^200 ■ We will examine here a possible deviation between electron and ion temperatures. These 
studies can help to understand how the cluster gas obtains hydrostatic equilibrium over large 
volumes. 

We define the average gas temperature as, 

(9) 

which will change over a typical electron-ion equilibration timescale, tci- We estimate the 
average gas temperature, kTg^s = Snl^^, by assuming a single power-law with 7 = 1.1 for 
the radial entropy profile, normalized in the cluster inner regions where Ti = because the 
relaxation times are much shorter there, figure 8(d) shows the ratio of the observed electron 
temperature to the estimated average gas temperature, where we have adopted Ui — 0.92 Ue 
(including He^+) for a fully ionized gas with X — 0.7 and Y — 0.28. Temperature inconsistency 
between Te and Tg^g is possibly larger than the simulation example (Rudd & Nagai 2009). 

The rapid Te decrease in the cluster outer regions may be explained by either the ICM not 
being in hydrostatic equilibrium or by differences between T^ and 7]. We could determine which 
interpretation is correct if we could directly estimate 7] from the line width. This measurement 
should be possible in the near future using the microcalorimeters on the ASTRO-H mission 
(Takahashi et al. 2008). 

5.5. Mass Estimation to r2oo 

We calculated the gravitational mass of A1413 to r2oo assuming spherical symmetry 
and hydrostatic equilibrium. From numerical simulations, these assumptions arc valid within 
~ 2 r2oo except for the core region at r < 0.3 r2oo, where cooling and heating of AGN give 
significant effects on the physical state of the gas (RoncareUi et al. 2006; Borgani et al. 2006). 
Previous X-ray studies mainly showed gravitational mass within rsoo because of instrumental 
hmitations. In this section, we determine the mass profile in the outer region of A1413. 

Assuming hydrostatic equilibrium, the total integrated gravitational mass, M<^r, within 
the 3-dimensional radius R is given by (Fabricant et al. 1980) 



kTR (dlnpg d\nT\ 



+ 



jirripG \d\n.R dluR 
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(11) 



where G is the gravitational constant, ji is the mean molecular weight of the gas and nip is 
the proton mass. We derive the above temperature and gas density profiles using the observed 
projected temperature and surface brightness profiles. We use the projected temperature di- 
rectly, but discuss the validity of this assumption below. We calculate the gas density from 
the normalization of the ICM spectral fit by taking into account the projection effect. The 
apec normalization parameter is defined as Norm = 10~^^/ ncn}idV/{4TT{l + z)'^D\) cm~^, with 
Da the angular diameter distance to the source. We estimated the de-projected nefin values 
assuming spherical symmetry and a constant temperature in each annular region as described 
in Appendix 1, and then assumed rie = 1.2 riu (excluding He^+) as described above. 

Allowing for the possibihty of Ti, we consider two cases for T: the electron tem- 
perature and the average gas temperature. We show the integrated mass profiles in figure 9(a) 
based on kT^ and kTgg^- These profiles are obtained without using any particular model since 
we perform the needed derivatives by differencing the temperatures and densities of adja- 
cent radial bins. The integrated mass within 13'.2'^q/^, which encompasses r2oo (14'. 8) is 
(8.8 ± 2.3) X IO^^Mq using kT^as- This mass is about 30% larger than that obtained using 
kTe of (6.6 ±2.3) X lO^^M©, although the difference is not statistically significant. The 30% 
difference in the temperatures propagates almost directly to the same mass difference. Our 
mass determination agrees with that of Vikhlinin et al. (2006), but not with Pointecouteau et 
al. (2005). These masses imply an overdensity with respect to critical of 177±47 and 132±47, 
where the errors are only from the mass errors. 

In the above mass estimation, we assumed that the observed projected temperature 
is the 3-dimensional value at the observed radius. Wc need to examine the systematic error 
caused by this assumption. In the following we denote the true 3-dimensional temperature of 
the ICM by T^^, which varies with radius. We derive the temperature from the spectral fit is a 
weighted mean of different temperatures projected along the line of sight. Often the projected 
temperature is defined as the emission- weighted temperature Te^, 



However, Mazzotta et al. (2004) discussed how the spectral response of an actual instrument 
implies that Tc„ can be quite different from what would be measured with that instrument 
observing a non-isothermal temperature distribution. For a better approximation, they intro- 
duced a spectroscopic-like temperature Tgi defined as. 



with a = 0.75, which empirically gave a good estimate of the T measured with XMM-Newton 
or Chandra. Rasia et al. (2005) and Shimizu et al. (2006) reported that the difference between 
Tew and Tgi can be as large as 30%. We carried out comparison of observed temperature with 
/cTew and kTgi in figure 9(a). The difference between kTe^ and kTgi takes the largest value of 



Jn^A{T)TdV 
Jn^A{T)dV ' 



(12) 



(13) 
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Fig. 9. (a) Comparison of observed temperature (black diamonds) with fcTew (doted diamonds), and 
kTsi (grey diamonds). Grey crosses show kT^ observed with XMM-Newton by Snowden et al. (2008). 

(b)Integrated mass profile (black diamonds: Suzaku with Tgas, grey diamonds: Suzaku with To, black 
crosses: XMM-Newton with Tgas, and grey crosses: XMM-Newton with Te). Vertical dashed line shows 
r200 = 14'.8. 

about 8.2% in the radius 2'.6 — 7'.0. These temperatures are consistent with the observed data 
with XMM-Newton. Taking a conservative value, our mass estimate would be more than 30% 
different from the true value because of our employment of the observed projected temperature 
as the 3-dimensional one. 



6. Summary 

• Northern outskirts of the relaxed cluster of galaxies A1413 was observed with Suzaku in 
the radial range of 2. '7 — 26' covering the virial radius of r2oo = 14'. 8. We excised 15 point 
sources above a flux of ~ 1 x 10"^"^ erg cm^^ s^^ (2-10 keV), and the CXB level after the 
point source excision was evaluated. We quantify all known systematic errors, and show 
statistical errors are dominant. 

• Suzaku detected X-ray emission of the ICM up to the 15' — 20' annulus beyond the virial 
radius. Significant temperature decrease to ~ 3 keV (factor of ~ 2) at r2oo is confirmed, 
which was reported in a few other clusters, PKS0745— 191 (George et al. 2008), A1795 
(Bautz et al. 2009), and A2204 (Reiprich et al. 2009). 

• Our entropy profile in the outer region (> 0.5 r2oo) joins smoothly onto that of XMM- 
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Newton at 0.15-0.5 r2oo, and shows a flatter slope of 7 = 0.90 ±0.12 than 7 = 1.1 (Voit 
2005) obtained with numerical simulations of adiabatic gas accretion. 

• Deviation of the entropy profile from the r~^^^ relation would show that electron temper- 
ature is not equal to gas temperature in outer region, where equilibration timescale for 
electron-ion collision, teu is longer than the elapsed time after the shock heating, ^elapsed- 

• The integrated mass of the cluster at the virial radius is approximately 7.5 x 10^^ Mq and 
varies by ~ 30% depending on temperatures (Te, Tgas, Te^, and Tgi) which we use. 

We are grateful to S. L. Snowden and A. Vikhlinin for communicating their unpub- 
hshed results. This work was supported by Grant-in-Aid for JSPS Fellows (20-2427) and the 

MEXT program "Support Program for Improving Graduate School Education" . JPH gratefully 
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